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Anisotropic flows (di, V2, us and Vi) of ligfit fragments up till the mass number 4 as a function 
of rapidity have been studied for 25 MeV/nucleon ''"Ca + "^^Ca at large impact parameters by 
Quantum Molecular Dynamics model. A phenomenological scaling behavior of rapidity dependent 
flow parameters u„ (n = 1, 2, 3 and 4) has been found as a function of mass number plus a constant 
term, which may arise from the interplay of collective and random motions. In addition, v/^jv-^ 
keeps almost independent of rapidity and remains a rough constant of 1/2 for all light fragments. 

PACS numbers: 25.75.Ld, 24.10.-i, 21.60.Ka 



Collective flow is an important observable in heavy ion 
collisions and it can bring some information on nuclear 
equation of state and in-medium nucleon-nucleon cross 
section [H, i, i i, i, i, 0, i, i, E [H El, [H, d. 
In intermediate energy heavy ion collisions, flows can 
reflect the interplay of collective and random motions. 
For a thermalized system, the random motions of emit- 
ted fragments are dictated by thermal energy, which is 
independent of mass. In contrast, fragment energy in- 
creases linearly with mass if there exists collective mo- 
tion [l^ [13, Ell • Many studies of the properties of the 
directed flow (v\) and the elliptic flow (v-i) have been car- 
ried out and much interesting physics has been demon- 
strated on the properties and origin of the collective mo- 
tion in both nucleonic or partonic levels. Recently, we 
carried out a Quantum Molecular Dynamics model cal- 
culation and found that there is a nucleon-number scaling 
for the transverse momentum dependent elliptic flow for 
light particles up to the mass number A = 4 [il|23,IIl|. 
But there are little studies about higher order flows such 
as U3 and U4 as well as their rapidity dependence. In this 
work, we shall present rapidity dependence of v\, V2, W3 
and W4 flows for a symmetric collision, namely ^''Ca -I- 
""^Ca collisions at 25 MeV/nucleon and large impact pa- 
rameters (6 = 4-6 fm) simulated by Quantum Molecular 
Dynamics model. The scaling behaviors of U4/(u2)^ as a 
function of rapidity is also presented. 

Anisotropic flow is defined as the different n-th har- 
monic coefficient u„ of the Fourier expansion for the par- 
ticle invariant azimuthal distribution 
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where (jj is the azimuthal angle between the transverse 
momentum of the particle and the reaction plane. Note 
that the z-axis is defined as the direction along the beam 



and the impact parameter axis which is labelled as x- 
axis. Anisotropic fiows generally depend on both particle 
transverse momentum and rapidity, and for a given ra- 
pidity the anisotropic flows at transverse momentum pt 

{pt = JpI+pI ) can be evaluated according to 



Vn[pt) = {cos{n(f)), 
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where (■ • ■ ) denotes average over the azimuthal distri- 
bution of particles with transverse momentum pt. The 
anisotropic flows Vn can further be expressed in term of 
single-particle averages: 
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where px and py are, respectively, projections of particle 
transverse momentum in and perpendicular to the reac- 
tion plane. 

The intermediate energy heavy-ion collision dynamics 
is complex since both mean fleld and nucleon-nucleon 
collisions play the competition roles. Furthermore, the 
isospin dependent role should be also incorporated for 
asymmetric reaction systems. Isospin dependent Quan- 
tum Molecular Dynamics model (IDQMD) has been af- 
filiated with isospin degrees of freedom with mean field 
and nucleon-nucleon collision [H, HI, [H, [H, [2||. The 
IDQMD model can explicitly represent the many body 
state of the system and principally contain correlation 
effects to all orders and all fluctuations, and can describe 



2 



the time evolution of the colhding system well. When 
the spatial distance (Ar) is closer than 3.5 fm and the 
momentum difference (Ap) is smaller than 300 MeV/c 
between two nucleons, two nucleons can coalesce into 
a cluster [2^. With this simple coalescence mechanism 
which has been extensively applied in transport theory, 
different size clusters can be recognized. 

In the model the nuclear mean-field potential is pa- 
rameterized as 
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where po is the normal nuclear matter density 
(0.16/m~'^), p„, pp and p are the neutron, proton and 
total densities, respectively, is z-th component of the 
isospin degree of freedom, which equals 1 or -1 for neu- 
trons or protons, respectively. The coefficients a, (3 and 7 
are parameters for nuclear equation of state. Cgym is the 
symmetry energy strength due to the density difference 
of neutrons and protons in nuclear medium, which is im- 
portant for asymmetry nuclear matter (here Csym = 32 
MeV is used), but it is trivial for the symmetric system 
studied in the present work. Vc is the Coulomb potential 
and C/^"'^ is Yukawa (surface) potential. In this work, 
we take a = 124 MeV, (3 = 70.5 MeV and 7 = 2 which 
corresponds to the so-called hard EOS with an inconi- 
pressibility oi K = 380 MeV. 

About 60000 events have been simulated for the col- 
lision system of ^°Ca -I- ^°Ca at 25 MeV/nucleon with 
impact parameter from 4 fm to 6 fm. The physics results 
were extracted at the time of 200 fm/c when the system 
has been in the freeze-out stage. Fig.[l](a) shows rapidity 
dependence of the average in-plane transverse momen- 
tum per nucleon (p^/A) around mid-rapidity for light 
fragments. Here rapidity y is the rapidity of fragment in 
the center of momentum (cm.) frame which is normal- 
ized to the initial projectile rapidity, i.e., y — ycmJUproj- 
It shows that (p^/A) decreases monotonously with the 
increasing rapidity, and decreases more rapidly for heav- 
ier fragments. The negative slope is mainly drive n by 
the attractive mean field in this energy region [ll|, |27| |. 
Fig. IHb) shows mass dependence of the so-called side- 
wards flow parameter (slope value) d{px/A)/dy, which 
can be easily extracted via linear fits near mid-rapidity 
in Fig. [Tl^a). The flow parameter is negative and its ab- 
solute value increases with mass apparently, which seems 
consistent with coalescence-like cluster production mech- 
anism 

In Fig. [21 we show integrated anisotropic flows wi , V2, 
V3 and V4 as a function of rapidity y. The trend of the 
directed flow vi shown in Fig. [UJa) is similar to that of 
(px/A), i.e., the directed flow decreases monotonously 
with the increasing rapidity around mid-rapidity. It is 
positive in target-like rapidity region and negative in 
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FIG. 1: (a) The average in-plane transverse momentum per nu- 
cleon as a function of the normalized rapidity in the c.m. frame, (b) 
Mass dependence of d{px/A) /dy. Squares represent for neutrons, 
circles for protons, triangles for fragments of A=2, diamonds for 
A=3 and stars for A=4. Error bars show the errors of the fits. 
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FIG. 2; (a), (b), (c), (d) show rapidity dependence of vi, V2, V3, 
V4 for light fragments respectively. Squares represent for neutrons, 
circles for protons, triangles for fragments of A=2, diamonds for 
A=3 and stars for A = 4. 



projectile-like rapidity region. And it gets to the ex- 
tremum near y = ±1.5, and then the absolute value of 
directed flow goes down with the increasing \ y \. It also 
shows at a given rapidity the absolute magnitude of vi 
is larger for heavier clusters. Fig. [Ifb) presents the ra- 
pidity dependence of elliptic flow V2- It shows that the 
flow is positive and the curve for a given fragment is sym- 
metric with y = 0, and descends monotonously when the 
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FIG. 3: (a), (b), (c) and (d) show rapidity dependence of di, V2, 
U3 and V4 per nucleon for light fragments, respectively. Squares 
represent for neutrons, circles for protons, triangles for fragments 
of A=2, diamonds for A=3 and stars for A=4. 



I y I increases. Again, heavier fragments have larger eUip- 
tic flow at the same rapidity. The rapidity dependences 
of vi and V2 are similar to the results at RHIC ener- 
gies for charged hadrons [2^ , but the interaction level of 
the matter in the two different energy region are obvi- 
ously different. In this low energy heavy ion collisions, 
the positive elliptic V2 essentially stems from the collec- 
tive rotational motion [O, H^]- Fig. [2Jc) and (d) show 
rapidity dependence of higher order flows of ^3 and V4, 
respectively. The trend of U3 and V4 are similar to that of 
vi and V2 , respectively, but they show smaller magnitude 
of flows for a given particle in comparison with vi and 
V2, respectively. Nevertheless the magnitude of flows of 
V3 and V4 are not too small and hence higher order flows 
should not be neglected in this energy region. 

For testing if the number of nucleon scaling works for 
rapidity dependent flows, we plot the flows per nucleon 
as a function of rapidity which is shown in Fig. [31 It 
seems that the curves for different fragments do not stay 
together as we imagined except for vi/A. One reason is 
that the effect of random motion which is independent 
of mass may weaken the mass scaling in rapidity space. 
For details, if we compare the magnitude of collective 
flows of other fragments with that of proton or neutron, 
it shows the flow magnitude of fragments with mass num- 
ber of 2, 3, and 4 is smaller than proton or neutron for 
vi and V2, but larger than proton or neutron for W3 and 
V4. However, if we use a mass number dependent func- 
tion C{A) = |(A-f |) instead of A for vi and V2, and use 
C{A) — 3(^— |) for V3 and V4, and plot Vn/C{A) versus y 
again which is shown in Fig.Hl we can now see the curves 
for different fragments coincide with each others excel- 
lently except for slight deviation in large rapidity region, 
i.e. in the spectators region. The coefficient C{A) may 
be qualitatively understood by two parts: the constant 
part reflects the effect of random motion which is inde- 
pendent of mass, and the part including mass number 



FIG. 4: Same as Fig. |3] but for v„/C{A) (n=l, 2, 3, 4) versus 
rapidity. 
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FIG. 5: The ratio of V4/V2 for neutrons (squares), proton (circles), 
the fragments of A = 2 (triangles), A = 3 (diamonds) and j4 = 4 
(stars) versus rapidity y. 



A reflects the collective motion which increases linearly 
with mass [l^, [2^ . From the mass dependent coefHcient 
of C{A), there are two classes: one for vi and V2, and an- 
other for V3 and W4. The larger coefficient for ^3 and 
may indicate that the contribution of collective motion 
to higher order flows W3 and W4 is larger than that to vi 
and V2. 

The RHIC experimental data demonstrated a scaling 
relationship between 2nd flow V2 and ?i-th flow w„ for 

hadrons, namely f„(pt) ~ v^.^'^iPt) aiidf4/i;| ~ 1.2 in the 
STAR data j3fl|] . And we have already shown the scaling 
of f4/w| as a function of pt in intermediate energy heavy 
ion collisions fH, H^, [2l| , i.e., the ratio oiv4,/v2 is inde- 
pendent of transverse momentum and shows a constant 
value of 0.5 for different light fragments. Similarly, we 
plot pt-integratcd v^/v^ as a function of rapidity which 
is shown in Fig. [5l The figure shows that the ratios of 
va/v2 for different fragments up to ^ = 4 roughly keep 
a constant of 0.5 in the studied rapidity region except 
for some ffuctuations at large rapidities. Therefore the 
ratio of W4/w| is almost independent of transverse mo- 
mentum and rapidity and its value is around 0.5. In the 
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frame of relativistic fluid dynamics which works in RHIC 
energies, the ratio value of 0.5 reflects that the collision 
system reaches to thermal equilibrium and its subsequent 
evolution follows the laws of ideal fluid dynamics [31] . A 
coincident value of 0.5 as relativistic fluid dynamics in 
the present intermediate energy may also indicate some 
kinds of thermal equilibrium has been also reached. 

To summarize, phenomenological behaviors of 
anisotropic flows as a function of rapidity for light 
fragments up till mass number 4 have been investigated 
by the simulation of 25 MeV/nucleon ^°Ga. + *°Ca 
collision in peripheral collisions in the framework of the 
Quantum Molecular Dynamics model. It was shown that 
vi and V3 of light fragments decrease monotonously with 
rapidity from positive value near target-like rapidity to 
negative value near projectile-like rapidity, while V2 and 
V4 are positive and show Gaussian-like shape with a 
peak around j/c.m. — 0. When we plot anisotropic flows 
per nucleon {vi/A, V2/A, V3/A and V4/A) versus rapidity 
for light fragments, the curves of different particles do 
not stay together perfectly except for vi/A. But when 
we plot Vn/C{A) (n = 1, 2, 3, 4) versus rapidity where 
C{A) is a linear function of the mass number A plus 



an additional constant term, the curves of different 
particles collapse onto the same curve. That indicates 
that the fragment flows may arise from the interplay 
of collective (a term proportional to mass number in 
C{A)) and random thermal motion (a constant term in 
C{A)) of nucleons. C{A) shows a classification between 
vi^2 and U3.4, this might reflect that the different role 
of collective motion on the different anisotropics in 
momentum space. Additionally, it was found that the 
ratio of v^/v^ are almost independent of rapidity and 
the value is about 0.5, which implies a thcrmalization 
scenario of the ideal fluid-like dynamics could be applied 
even in such low energy heavy ion collisions. Further 
theoretical investigations and experimental checks are 
awaiting. 
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